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Abstract 

A simple approach is proposed for the quantization of the electromag- 
netic field in nonlinear and inhomogeneous media. Given the dielectric 
function and nonlinear susceptibilities, the Hamiltonian of the electromag- 
netic field is determined completely by this quantization method. From 
Heisenberg's equations we derive Maxwell's equations for the field opera- 
tors. When the nonlinearity goes to zero, this quantization method returns 
to the generalized canonical quantization procedure for linear inhomoge- 
neous media [Phys. Rev. A, 43, 467, 1991]. The explicit Hamiltonians for 
the second-order and third-order nonlinear quasi-steady-state processes are 
obtained based on this quantization procedure. 

PACS numbers:42.50.-p, 42.65.-k, ll.lO.Lm 



*E-mail:gcguo@sunlx06. nsc.ustc.edu.cn 



1 



1 Introduction 



Early quantization of the electromagnetic field is performed in empty cavities or 
in infinite free space [1] . However, with the growth of interest in quantum optical 
phenomena taking place inside material media, several approaches have been pro- 
posed for quantization of the electromagnetic field in nonlinear, inhomogeneous, 
or dispersive media [2-21]. Early attempts towards quantization of the nonlin- 
ear media, while incorporating the known linear theory, did not fully reproduce 
the nonlinear field equations [4]. An innovative treatment was first proposed by 
Hillery and Mlodinow who successfully quantized a nonlinear medium by intro- 
ducing the dual potential [6]. Later, Drummond extended the Hillery- Mlodinow 
procedure to dispersive media [9] . There are also other approaches in this direc- 
tion. Glauber and Lewenstein generalized the canonical quantization method by 
modifying the gauge condition to deal with the inhomogeneous linear media [11]. 
Abram and Cohen, following the canonical quantization procedure, presented a 
quantum formulation for light propagation in nonlinear effective media [12]. And 
recently, Santos and Loudon gave an alternative approach to the quantization of 
the electromagnetic field in linear one-dimensional dispersive media [19]. Devel- 
opments towards the absorbing dielectrics also appeared [20,21]. 

In this paper, we propose a relatively simpler approach to the electromagnetic 
field quantization in nonlinear and inhomogeneous media. The procedure follows 
Ref. [11] in using the material independence of the commutation relations for 
the fields D and B , pointed out by Born and Infeld [2], as a starting point in 
the quantization. We extend this to nonlinear media in which D and B can 
be expressed as isochronous functionals of the fields E and H . D and B are 
expanded into the mode functions. Furthermore, we explicitly derive Maxwell's 
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equations for the field operators from Heisenberg's equations. This procedure is 
apphed to the quantization of the second-order and third-order nonhnear quasi- 
steady-state processes and we obtain the exphcit Hamiltonians. Though these 
Hamiltonians are aheady in wide use in quantum optics, their derivations are 
mainly based on the early quantization procedure by Shen [4,22,23] and known 
by now to be inconsistent with the nonlinear field equations [9] . So here we give 
a justification of these Hamiltonians. 

The arrangement of the paper is as follows. The quantization procedure is 
proposed in Sec.2. Given dielectric tensor and nonlinear susceptibilities, this 
quantization procedure completely determines the Hamiltonian of the electro- 
magnetic field, which is expressed by annihilation and creation operators. Then 
we derive Maxwell's equations from Heisenberg's equations for the field operators. 
In Sec. 3 we show this quantization procedure returns to the generalized canon- 
ical quantization method in Ref. [11] when the medium is linear. The explicit 
Hamiltonians of the second-order and third-order nonlinear quasi-steady-state 
processes are obtained in Sec. 4 by application of this quantization procedure. 

2 Quantization in the presence of nonlinear me- 
dia 

We consider the electromagnetic field in nonlinear media, which may be inhomo- 
geneous. The source-free Maxwell equations in matter take the forms [24] 



V • D = 0, 



(1) 



V • 5 = 0, 



(2) 



IdD 



V X H, 



(3) 



c dt 
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In nonlinear media, D (t) and B (t) are complicated nonlinear functionals of 
E (t) and H (t). From the Maxwell equations, the energy density U of the 
electromagnetic field in nonlinear media is determined by 

dU (t", = E (t", t) ■ dD (l^, t) + H (l^, t) ■ d'B {j, t) . (5) 

The Hamiltonian (or the energy) is 

H = j d^T'uij.t) . (6) 

For the electromagnetic field in linear media, the canonical quantization method 
is generally used. The vector potential A is chosen as the general coordinate and 
the Columb gauge V ■ A = is often used. A and its conjugate momentum can 
be expanded into a set of transverse complete spatial functions and the expansion 
coefficients are expressed by annihilation and creation operators. Then substitut- 
ing the expansions of A and its conjugate momentum into the Hamiltonian, one 
achieves quantization of the electromagnetic field in linear media. However, for 
the nonlinear media, the canonical quantization becomes much more involved. 
A and E were chosen as the canonical variables in the early treatments, which 
did not incorporate Eq.(l). In fact, no rigorous approach had been proposed for 
nonlinear media until Hillery and Mlodinow introduced the dual potential and 
then followed the canonical quantization procedure. Here, inspired by the result 
in Ref. [1] that the fields D and B have medium-independent commutation rela- 
tions in inhomogeneous linear media, we choose the fields D and B , rather than 
or A , as the starting point of the electromagnetic field quantization. This 
choice is also consistent with the results in Ref. [6] [9] and [12], where the field D 
was found to be the canonical momentum. Starting from the mode expansions 



of the fields D and B , we can present a concise formulation of the quantization 
and a clear derivation of Maxwell's equations for the field operators. 

Prom Equations (1) and (2), the fields D and B can be expanded into a set 
of transverse complete spatial functions | / (^)} ^^"^ ^ T (^)} ' 
respectively, 

s(^.*)=-i:^T.w7v„(^). p) 



k /i 



B {l^,t 



k n 



The expansion functions and coefficients satisfy Hermitian conditions 



(8) 



(9) 
(10) 

(11) 

In addition, the functions / {/^^ satisfy transversality, orthonormality and 
completeness conditions 



f * 
k // 


= / 


— k fi 


^ k fx 






k fi 


= p 


k /i 



i:7tJ^)7t„,(^')=4(^-^'). 



(12) 
(13) 
(14) 



k IX 



where the transverse (5— function 5t is defined as 



1 



(T k 



Si 



^3 



\ 



i k ■ r 



(15) 



The transversality condition (12) makes the completeness equation of / 



k /i 



take the form of (14). In free space, the plane wave is chosen as the expansion 



function 



where the unit vectors (a* = 1)2) satisfy 

The expansions (7) and (8) have the same forms as those in hnear media. 
We further suppose that the expansion coefficient operators satisfy the same 
commutation relations. So 



ihS-.-.,5,„/. (18) 

k k' l^i^ ^ ' 



The fields E {jr* , and H ^7^, can be expressed by D (^'r^, and B ^T^, 
from the nonlinear functional relations between them. Prom (5) and (6) the 
Hamiltonian of the electromagnetic field becomes a nonlinear functional of D (t^, tj 
and B (T^,tj. After substituting the expansions (7) and (8) into it , we get the 
Hamiltonian, which is expressed by annihilation and creation operators. Given 
the functional relations between E (l^,tj,H (^5^) and D (l^,tj,B (y^,tj, 
the Hamiltonian form is completely determined by this quantization procedure. 

Now we show reasonableness of the quantization method. Comparing (7), (8) 
and (18) with the corresponding equations in Ref.[l], we know that when the 
nonlinearity goes to zero the above procedure returns to the generalized canonical 
quantization method. Furthermore, this quantization gives the correct Maxwell 
equations. In the following we derive Maxwell's equations for the field operators 
from Heisenberg's equations. 

Prom the transversality of the expansion functions / , the first two Maxwell 

k fi 

equations (1) and (2) are obviously satisfied. Equations (7) (8) and (14) (18) give 
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the commutator of the field operators D and B 



Di (r , t) , Bj (r', t)] = the (V X 5) . . ( 



(19) 



In the derivation, the following relations are used. 



X ^")^, (^) =Sn '^mndr^J^ J d^t (^6., - 



\ k -T* 



— Z-( '-■imn'-'m (27r)'^ "^nj"^ 



(V X 5),^. (r) 



and 



(20) 



(V' X 5),, (j - T*') = - (V X 5) .. (t* - T") 

mn ^ ' 

= - E SjmidmS (l^ - 1^') (21) 

=E £imjdmS (T^ - T"') 

= (V X 5)^^. (r - r') . 

The commutator (19) has been given in Ref. [12] in one- dimensional form. Here 
we extend it to the general case. From the commutator (19), the commutator 
between D or B and an arbitrary functional F (^D , B^ of D and B , which 
may be nonlinear, can be expressed by functional derivation as follows: 



Di(T',tj ,F = ihc Eimndr, 
mn 



5Bn(j,t) 



-F. 



(22) 



(23) 



Prom these two equations and (5), (6), the Heisenberg equations of the field op- 
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erators D and B take the forms 



1 dPj _ J_ 
c dt ihc 



mn SBn[ r ,t\ 



-L^SimnOmg^^^^^^ (24) 

Z_(Z_( '-'imn'-'m-'^-'^j f)R 
mn j " 



V X H 



-iWx E] . (25) 



Similarly 

c dt 

So we have clearly derived Maxwell's equations for the field operators from Heisen- 
berg's equations. The isochronous commutators (22), (23) play an essential role 
in the derivation. The derivation holds for linear or nonlinear media. However, 
for dispersive media, nonlocal relations in time between the Hamiltonian H and 
the fields D , B arise and the isochronous commutators (22), (23) cannot be ap- 
plied in this case. Here we meet the long standing difficulty in quantum optics 
in quantizing nonlinear and dispersive dielectrics. 

3 Quantization of the electromagnetic field in 
linear inhomogeneous media 

In this section we use the above method to quantize the electromagnetic field in 
linear inhomogeneous media. The medium is characterized by 

A (t^, t) e,, (r) (r , t) , (26) 
j 

'B {r,t) = H {r,t). (27) 



The Hamiltonian (or the energy) is 

(28) 



In free space the expansion function f is expressed by Eq.(16). Substi- 
tuting the expansions (7) and (8) into Eq.(28) and noting 

k fx , 

'huj-. 



where cu-^ = k 

k fj, 

creation operators. 



c , we get the Hamiltonian expressed by annihilation and 



X 



(31) 



V* , (l^, t') at at, -V*,( t, -~k') at a^, , + h.c. 
where V„,i is defined by 

If the second-order tensor e"^ is a scalar, "e*^ ■ ~e ^, , in Eq.(32) can be put 

k fi k fi i\/ 

out of the integration. This case had been discussed in detail in Ref.[ll] by the 
generalized canonical quantization method. Here we get the same results. 

4 Quantization of the quasi-steady-state optical 
field in nonlinear media 

In this section we consider quantization of the quasi-steady-state optical field 
in nonlinear media. The second-order or third-order nonlinear process is most 
important [25,23]. First we apply the quantization method to the parametric 
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process. The optical field is composed of three quasi-monochromatic fields with 

central frequency ui,U2,ui3, respectively, and LO3 = ui + U2, i.e., 



'E (t) =^ 'E^^ (t) e~^'* + h.c, 



(33) 



where E (t) is slowly varying amphtude. Under the quasi-steady-state approx- 
imation, the term E^^'^ (t) e"*"^'* in Eq.(33) can be viewed as a monochromatic 
field with frequency cui, i.e., the dispersion of the optical field in the medium is 
negligible. Suppose the refractive index is independent of orientation of space. 
Then the field D^^^ (t) can be expressed by E^^^ (t). For example, 

D (t) = (us) E^'^ (t) + x^') (cos = a;i + 002) : 'E^'^ (t) 'E^^^ (t) , (34) 

where (cos) is introduced phenomenologically to show the component of the 
optical field E (t) has been viewed as a monochromatic field with frequency 
ujs- D {t) and D (t) have similar expressions. From these equations, E [t) 
(i = 1, 2, 3) may be expressed by D [t] as 

E it) = ^7-^ - 7^'^ (^3 ^u,-rU2):D (t) D (t) , et.al., (35) 

n [UJ3) 



where 



1 



(tjg ^uJi+ U2) 



(36) 



In the derivation of Eq.(35) the approximation that the nonlinear terms are much 
smaller than the linear terms is used. Substituting (35) into (5) and (6) and 
making the rotating wave approximation, we get the Hamiltonian of the electro- 
magnetic field, which is expressed by D and B 



3 

E 



^ + Lb « it) e-^^^ + h.c 



(37) 



-2 



7(2) (a;3 = cui + UJ2) ;^(')+ (t) ^(2) (t) (t) + h.c. | 
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In the derivation, the holo-exchange symmetry of the tensor x^'^^ has been used. 

In the expansions (7) and (8) of the fields D , B , only the terms with the 



subscripts 



k 



ki, where ki 



n{uJi)uJi 



,{i — 1,2,3), make contributions to the 



interaction. So the expansions can be simplified to 



D « (t) : 



—iuiit 



= ^1 



y yy^a^ (t) , 



(38) 



B (t) e"^'"'* = c\ 



n 



(39) 



The function / -> can be decomposed as f ^ = ( r 



f-g (^'r*^ satisfies the eigen-equation 



k 



e -+ , where 



(40) 



Substituting (38) and (39) into Eq.(37), we get the Hamiltonian expressed by 
annihilation and creation operators 



ki =ki 
k i,k2,k 3 



E 



J2 fioJiO^ a—> 

i=l k iHi k iHi 



+ («-» -+ -+ ai^ a-* a-+ + h.c.) 

k iixi k 2At2 k 3/i3 k 1 k 2 k 3 fc g^^g fe 2/U2 iMi 

where the constant (3 (called phase-matching factor) is defined as 



(41) 



13- 



k 1 k 2 k : 



J k 3 k 2 k 1 



(42) 



and V is volume of the nonlinear media. The constant a is 



a- 



k lfJ.1 k 2M2 k 3^13 



/ h^^^^u)2^^3 



8V')i2(Lji)n2(a)2)n2(u;3) 



fc*3M3 2M2 k 



(43) 
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and it is determined by the polarization matching condition of the optical field. 
When k I, k 2, k ^ are coUinear, the Hamiltonian (41) can be simplified. Sup- 
pose the polarizations of the optical fields are given (indicated by fj!^, fj,^, fJ-'s 
spectively) and k j = ki~e^z (^ — 1)2,3), where ~e'z is the unit vector of z-axis. 
The expansion function is approximately a plane wave, i.e., 

f,^(r)=S,^{x,y)^, (44) 

where L is the interaction length. Let represent the operator a, / and 

" = ^k^^^[k,^^',k3^.',\[J j Sl^Sk,Sk,dxdy, (45) 
then the Hamiltonian (41) is simplified to 

H =Y^ fiUiajai + \ a ■ ^^^^ 0^0201 + /i.e. j , (46) 

where the phase mismatch A A; = ki + k2 — k^. Eq.(46) is often used to analyze 
quantum properties of the parametric process [26-28]. Here we give its exact 
derivation and determine the expression of the parameter a. 

The optical field in the third-order nonlinear medium can be quantized in 
a similar way. For example, the Hamiltonian of the nondegenerate four-wave- 
mixing process [23] with ujz -\- oJi — oJi + 0J2 is 



k i — k-i r 4 

V V V /^l'M2,/i3,/i4 
Kl,fc2,'C3,fc4 



E fiUioL, a-^ 

.1=1 k iixi kiHi 

(47) 



/^T* -r -r «if '^if + ^-c-) 



where 



^fcifc2fe3fc4 J ki ks'' k2-' ki ' ^ ' 
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and 

^ ^ = —i / n*uiiu!2U!3W4 

"fc i/ii'fc2/i2'fc3/i3'fc4Ai4 V 16V^n^(a;i)n^(a;2)n^(w3)n'-^(a;4) 

(49) 



k AliA k 3/U3 * 2^*2 fc l/il 

The Hamiltonians (41) and (47) make the foundation for analyzing quantum 
properties of the parametric or four-wave-mixing process. 
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